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ABSTRACT 

In this paper, we report a refined determination of the orbital parameters and the detection of the Rossiter-McLaughlin effect of the 
recently discovered transiting exoplanet HD147506b (HAT-P-2b). The large orbital eccentricity at the short orbital period of this 
exoplanet is unexpected and is distinguishing from other known transiting exoplanets. We performed high-precision radial velocity 
spectroscopic observations of HD 147506 (HAT-P-2) with the new spectrograph SOPHIE, mounted on the 1.93 m telescope at the 
Haute-Provence observatory (OHP). We obtained 63 new measurements, including 35 on May 14 and 20 on June 11, when the planet 
was transiting its parent star. The radial velocity (RV) anomaly observed illustrates that HAT-P-2b orbital motion is set in the same 
direction as its parent star spin. The sky-projected angle between the normal of the orbital plane and the stellar spin axis, A = O.l^'i^ l °, 
is consistent with zero. The planetary and stellar radii were re-determined, yielding Rp = 0.951^|] J[jj Rj„p, Rs = 1.416^J[|]*2 Rq. The 
mass (Mp = 8.62^j[ jj Mj„p) and radius of HAT-P-2b indicate a density of 12. 5^^^ g cm"', suggesting an object in between the known 
close-in planets with typical density of the order of 1 g cm"', and the very low-mass stars, with density greater than 50 g cm"\ 
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1. Introduction termediate object between giant planets and low-mass stars, near 

IJ ■ the brown dwarf population. Its der isity is indeed close to the up- 

. ^ ; Thirty of the almost 270 known exttasolar planets have been pgr limit of the planetary models (iBaraffe et al.ll2003l) and may 

k> observed to transit their parent star^ This small subgroup of put this object in the transition region between planets and brown 

5—1 planets have the highest impact on our overall understanding of dwarfs 

C3 close-in giant planets because we can estimate their mass and We report here new RV measurements of HD 147506 (HAT- 

radius, and for some of them directly study their atmosphere, p.2) obtained to provide additional information and constraints 

These transiting hot Jupiters have masses from 0.07 to about 8 on this unusual planetary system. This was done by 1) refin- 

Mjup and radii from 0.3 to about 1.4 Rj,,,. This set of planets ing the orbital parameters, 2) refining the mass and radius of 

was recendy extended to hot Neptunejjbjects by the discov- the companion and 3) measuring and modeling the Rossiter- 

ery of the transit of GJ436b (|Gillon et al., 2007), and to super- McLaughlin effect (RV anomaly due to the partial eclipse of the 

massive planets by the detection of HD 147506b (HAT-P-2b) by rotating stellar surface). We present here a more precise orbital 

iBakos et al.| WB) (hereafter B07). The discovery paper of the solution and the measurement of the sky projection of the incli- 

latter derived the key orbital and physical parameters of this ex- nation between the stellar spin axis and the normal of the orbital 

ceptional object, which differs by its mass (9.04 My„p), orbital plane 
period (5.63 days), and eccentricity (0.52) from other transit- 
ing planets. The radius of HD 147506b implies an uncommon 
measured mean density (11.9 gem"-') and surface gravity (227 

ms"'). These properties suggest that HAT-P-2b might be an in- The parent star of HAT-P-2b, HD 147506, was observed in 

May and June 2007 with the new spectrograph SOPHIE 

Send offprint requests to: B. Loeillet, e-mail: benoit.loeillet@oamp.fr dBouchy & the Sophie tea^|2006|) mounted on the 1 .93-m tele- 

* Based on observations collected with the SOPHIE spectrograph on scope of Haute Provence Observatory. Observations were con- 

the 1. 93-m telescope at OHP, France (programs 07A.PNP.MAZE and ducted in the high-efficiency mode (HE mode), which provides 

07A.PNRCONS) a spectral resolution of R~ 39,000. The Thorium-Argon lamp 

' see .http://obswww.unige.ch/~pont/TRANSITS.htm| was used to calibrate the wavelength scale. The simultaneous 
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ThAr calibration mode was not used, as wavelength calibration 
was performed less than 1 hour before and after the observations 
leading to an instrumental stability of less than 2ms"'. 

We obtained eight out-of-transit spectra on six nights be- 
tween May 6 and 15, 2007, and a sequence of 6 hours, com- 
prised of 35 spectra during and after the transit on May 14 
(JD=2454235). Unfortunately, we could not observe the ingress 
phase of the transit that occurred just before twilight. Another 
sequence was obtained during transit on June 1 1 (JD=2454263) 
to get a full coverage of the transit. The typical exposure time 
was 10 minutes, long enough to reach a sufficiently high signal- 
to-noise ratio (SNR) and short enough to adequately sample the 
observed transit. The average SNR per pixel, at A =5500 A, was 
about 70 during the first observed transit and about 90 during the 
second one. Another limitation of high-precision RV measure- 
ment in the mode used (HE) concerns the guidance centring of 
the target in the fiber. As the diameter of the fiber is quite large (3 
arcsec), a strong decentring under very good seeing conditions 
(< 1 .5 arcsec) could indeed induce a RV shift of a few tens m s"' . 

We determined RVs using a weighted cross-correlation 
method, followin g the procedure of Baranne et al. ( 199^ and 
iPepe et all (l2005l) . with a numerical mask constructed from the 
solar spectrum atlas corresponding to a G2 dwarf star. This stan- 
dard mask is well adapted to the F8-type spectrum of the primary 
star. We estimated the measurement uncertainties based on the 
pho ton-noise empir i cal rel ation detailed by Bouchv et al. (l2005|) 
and ICameron et all (l2007l) . For the spectra obtained during the 
night of the first and the second transit, the typical photon-noise 
uncertainty is 14 and 10 ms"', respectively, whereas for the 
other 8 spectra this uncertainty ranges from 6 to 28 ms"'. The 
journal of SOPHIE observations, including barycentric Julian 
dates (BJD), RVs, photon-noise uncertainties, and SNR per pixel 
is given in Table [T] (online material). The phase-folded RVs are 
plotted in Fig. [1] As illustrated by the higher SNR values for 
the data acquired during the second night of observations, we 
have strong evidence that the observational conditions were very 
good. As explained above, the data are thus aff'ected by guiding 
noise and we assumed an additional systematic error of 50 m s"' . 
We estimated this systematic error thanks to later instrumental 
tests under the same conditions. 

3. Stellar properties 

As described in B07, the spectroscopic determination of the ra- 
dius of HD 147506 is very sensitive to the method used, as well 
as to the log g determination, which may be aff'ected by the un- 
certainty of the spectrum continuum and the large projected ro- 
tational velocity, v sin I, of the star. The spectroscopic approach 
described in B07 provides a stellar radius of 1.474^y '|g^ R©. Our 
combined SOPHIE spectrum provides an independent spectro- 
scopic determination of the stellar parameters (Rs - 1-416^q q^2 
Ro), consistent with B07. 

^From the full width half maximum (FWHM) of the av- 
eraged cross-coiTelation functions (CCF) of SOP HIE spectra, 
which were calibrated to yield stellar v sin /, values dSantos et al.l 
12003 ). we determined the vsin/, of HD 147506 to be equal 
to 21.3+1.3 kms"', somewhat slightly larger than the veloc- 
ity derived by B07 (19.8+1.6 kms"'). The metallicity index 
[Fe/H1=0.l 1 +0.10 we obtained, using the method described in 
ISantos et al.l (|2002|) is in full agreement with B07. The activ- 
ity index of log R'^ ^--4. 15 +0.02 was derived from the H and 
K Call lines and appears to be close to the value determ i ned b y 
B07 (log/?;y^=-4.72+0.05). As estimated in lSantos etaTlfcOOOh . 
our activity index for a F8 dwarf star implies a stellar jitter from 



Table 1. SOPHIE Radial velocities of HD147506 (HAT-P-2) 
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0.0148 


66 


54231.5987 


-19.5311 


0.0121 


78 


54235.3466 


-20.1916 


0.0156 


61 


54235.3538 
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0.0180 
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-20.3008 


0.0167 
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0.0209 


46 


54235.3938 
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0.0163 
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-20.4530 


0.0180 


53 
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0.0150 
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0.0144 
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-20.5032 


0.0157 


60 


54235.4456 


-20.5318 


0.0178 


53 


54235.4535 


-20.5340 


0.0182 


52 


54235.4608 


-20.4990 


0.0143 


66 


54235.4681 


-20.5047 


0.0124 


75 


54235.4759 


-20.4747 


0.0117 


81 


54235.4831 


-20.4974 


0.0113 


82 


54235.4904 


-20.5364 


0.0111 


84 


54235.4981 


-20.5384 


0.0111 


85 


54235.5054 


-20.5484 


0.0108 


87 


54235.5126 


-20.5785 


0.0135 


70 


54235.5204 


-20.6029 


0.0142 


66 


54235.5277 


-20.6066 


0.0153 


62 


54235.5350 


-20.5946 


0.0115 


81 


54235.5434 


-20.6155 


0.0118 


80 


54235.5507 


-20.5987 


0.0113 


82 


54235.5580 


-20.6266 


0.0110 


85 


54235.5682 


-20.6683 


0.0119 


80 


54235.5755 


-20.6712 


0.0123 


77 


54235.5827 


-20.6585 


0.0115 


82 


54235.5905 


-20.6819 


0.0108 


87 


54235.5978 


-20.7002 


0.0128 


74 


54235.6051 


-20.7205 


0.0117 


81 


54236.5190 


-20.2207 


0.0056 


99 


54263.4521 


-20.2050 


0.0157 


71 


54263.4594 


-20.1880 


0.0100 


102 


54263.4666 


-20.2048 


0.0089 


113 


54263.4739 


-20.1887 


0.0092 


110 


54263.4804 


-20.1626 


0.0114 


90 


54263.4860 


-20.1750 


0.0120 


86 


54263.4915 


-20.1832 


0.0128 


81 


54263.4971 


-20.1429 


0.0096 


104 


54263.5030 


-20.1565 


0.0102 


99 


54263.5086 


-20.1615 


0.0108 


95 


54263.5141 


-20.1610 


0.0104 


97 


54263.5236 


-20.1810 


0.0103 


99 


54263.5291 


-20.1955 


0.0151 


72 


54263.5347 


-20.3046 


0.0179 


63 


54263.5411 


-20.3391 


0.0277 


48 


54263.5561 


-20.3457 


0.0127 


85 


54263.5634 


-20.4009 


0.0133 


82 


54263.5706 


-20.4169 


0.0173 


67 


54263.5779 


-20.4377 


0.0107 


99 


54263.5852 


-20.4763 


0.0091 


113 
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a few ms 'to ab out 20 ms~', w hich is confirmed by the cal- 
ibration made by 'Wrighj (|2005|) . Analysis of the line-bisector 
computed for all out-of-transit spectra does not show significant 
variations nor correlation with the RVs. 



4. Determination of the planetary system 
parameters 

In this section, we first describe the procedure used to fit the 
observed velocities with the Keplerian orbit and the RM RV 
anomaly, and how we estimated the uncertainties. In the second 
and third subsections, we discuss our results. 



4.1. Analysis of the RV data 

We used all available high-precision spectroscopic data to model 
the orbit and the RM effect simultaneously. This data includes 
10 Lick spectra and 13 Keck spectra obtained by B07, and the 
63 SOPHIE spectra. Our model is comprised of 15 parameters: 
The period, P; periastron passage time, Tq; orbital eccentricity, 
e; angle between ascending node and periastron, co; RV semi- 
amplitude, K; RV zero point, Vo (these first six are the classi- 
cal orbital parameters); planetary to stellar radii ratio, Rp/Rs', or- 
bital semi-major axis to stellar radius ratio, a/Rs', angle between 
sky projection of the orbital angular momentum axis and stel- 
lar spin axis. A; line of sight stellar rotational velocity, vsin/,; 
orbital inclination angle, /; and the stellar linear limb darken- 
ing coefficient, e. We have also determined a velocity shift be- 
tween Keck and Lick velocity zero points, Avkl, and SOPHIE 
and Keck zero points, Avsk- We also estimated a velocity shift, 
Avs2, between the SOPHIE measurements taken on the sec- 
ond transit night (June 11) and the rest of the SOPHIE mea- 
surements, taken about a month earlier The period was fixed 
on the value given by B07 (P - 5.63341 days), consider- 
ing its very high accuracy (lis) derived from extensive photo- 
metric observations. A linear limb-darkening coefficient e=0.71 
was also used, considering the stellar T^ff — 6250 K jClaretl 
12004 . Hence, our model has 13 free parameters, where 8 are 
non-linear (To, e, a>, Rp/Rs, a/Rs, A, v sin I,, and /) and 5 linear 
(K, Vo, Avkl, Avsk and Avsi)- 

We used a Keple rian model for the orbit, and the analytic ap- 
proach described b ylOhta et al.l (l2005i) . to model the RM effect. 
Equations given by lOhta et alT ( l2005l) for the RM RV anomaly 
were modified to make them dependent on Rp/Rs and a/Rs, in- 
stead of Rn,Rs and a. 



HD147506 (HAT 



SOPHIE 



We searched the parameter space for the global minimum 
position, using the equation below. We used the linear least 
square method, along with the data, to calculate the linear pa- 
rameters at each position in the parameter space of the non-linear 
parameters. We modified omx^ function to account for external 
information, namely the line-of-sight stellar rotational velocity, 
derived here from the spectra, the radii ratio, and transit duration, 
fromB07: 



RVo,i-RVcj^2 ,vsin/,-21.3 , 

-) +( 7^^ ) 



RV, 



Rp/Rs-QMU 
+( ^^^^^ ) +(- 



1.3 



•0.177, 



0.0009 0.002 

where RVo.i and RVc,i are the /-th observed and calculated RVs 
and RVfrr,i is its error. T^,,,. is the transit duration and is related 
to a/Rs, Rp/Rs, i and the orbital parameters P, e, and oj. The un- 
certainties were computed directly from the linear least squares 




„ 0.3 

"„ 0.2 

S °l 

^ 



-0.2 
-0.3 



, 1 , 1 , 1 , 1 , 1 , 1 , 1 , T 1 " 


0.1 


0.2 0.3 0.4 0.5 0.6 0.7 0.8 


0.9 1 


- 1 ' 

- 1 , 


1 1 1 1 1 i_ 1 1 1 1 1 1 




1 ' 1 - 

1 V 1 - 



0.1 O.a 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Orbital 



Fig. 1. Top: Phase-folded Radial Velocity measurements of 
HAT-P-2 superimposed on the refined Keplerian orbital solution. 
Open circles refer to the Keck and Lick measurements. Filled 
circles refer to the SOPHIE measurements. The inset shows 
a zoom around the transit where radial velocities exhibit the 
Rossiter-McLaughlin effect. Bottom: Residuals around the or- 
bital solution. 



analysis for the linear parameters ("Press e t aLlll992l Sect 15.4, 
Eq 15.4.15). The procedure is a bit more complex concerning 
the non-linear parameters. For each specific parameter, we in- 
crease and decrease it using small steps starting from the mini- 
mum solution value. At each step, the rest of the parameters 
are fitted while holding the specific parameter constant. Then, a 
4th degree polynomial is fitted to the x^ values obtained and the 
1-cr uncertainty is estimated by identifying the values of the fit 
corresponding to the minimum value -i-l. 

In the fit procedure, we adopted the B07 stellar jitter of 60 
ms ' for the Lick and Keck measurements. For the SOPHIE 
measurements, we adopted a jitter of 17 ms"', resulting in a re- 
duced = 1 for the out-of-transit SOPHIE measurements. This 
value is in a good agreement with the value estimated from our 
revised activity index and should be compared to the jitter of 60 
ms ' mentioned in B07, considering the span covered by Lick 
and Keck observations was 240 nights. Over such a long span, 
one cannot exclude variations in the stellar activity level. We 
also note that the effect of this jitter may have a time scale com- 
parable or longer than the time of a single exposure, inducing a 
correlated noise effect. 



4.2. Results of the fitted orbital solution 

Table |2] lists the result for the fitted parameters of our model. 
The solution depends on the RV data and the prior constraints 
(vsinls, Rp/Rs, transit duration), derived from the spectroscopic 
analysis in this work and the photometry from B07. The re- 
fined orbital solution is plotted in the upper panel of Fig.[T] and 
Fig. |2]presents the RM RV anomaly model after subtraction of 
the Keplerian orbit. 

The RMS of the out-of-transit residuals of SOPHIE mea- 
surements, spanning 36 days, is equal to 18 ms"', whereas 



4 



Loeillet B. et al.: Refined parameters and spectroscopic transit of HD147506b 



RMS of all out-of-transit measurements, spanning 282 days, is 
50 ms"'. The orbital parameters we derived are consistent with 
B07. However our uncertainties are smaller, due to a larger sam- 
ple of high-precision RV data. 

We searched for a second planetary signal in the RV residuals 
to look for hints of a third body in the system. No clear periodic 
signal appears in the RV residuals. We estimated that we can 
exclude the presence of a second planet of mass greater than 1.3, 
1.5 and 1.8My„p for an orbital period shorter than 50, 100, and 
200 days respectively. However, the increased RMS for all out- 
of-transit residuals suggests that a long-term RV follow-up of 
this star is needed. 



4.3. Measurement of the Rossiter-McLaughlin effect 

As an important result, the sign of the RV anomaly shows 
that the orbital motion is set in the same direction as the stel- 
lar spin, similar to the four previous observed RM effects on 
transiting exopl anets (O uelozetal. (2000); Winn et al. (2005b 
for HD 209458; IWinn et alJ (20061) for HD18 9733; IWolf etaLl 
(120071) for HP 149026: and lNarita et al] (|2007|) for TrES-1). The 
orbital inclination angle / we derived is in full agreement with 
the value determined by B07. The projected rotation velocity of 
the star vsin/, determined from the RM fit is 22.9+1.2 kms"'. 
This value is greater than our spectroscopic determination from 
SOPHIE CCFs, and 2-cr gre ater than the determination of B07. 
However. IWinn et al.l (|2005 |) showed that v sin 7, measured with 
the analytical formulae from Oh ta et alj (|2005) is biased toward 
larger velocities by approximately 10%. Moreover, as the planet 
crosses the star at the equatorial plane, assuming a differential 
rotational velocit y of the star , the fitted v sin /, coiTesponds to the 
maximum value. [Ktiker & R iidiger (2005) show that a differen- 
tial rotational rate can be as high as about 10% for a F8-type star. 
The projected angle between the stellar spin axis and the normal 
of the orbital plane, A - 0.2 ± 12.5°, is consistent with complete 
alignment of the stellar spin and orbital angular momentum. 

We extended our analysis in the search for the global min- 
imum without any constraint on the vsin/^. The result ob- 
tained gives a significantly larger v sin /, value (29.5 '^22 
and a consistent lambda value (5.0^^^^ °). However, such a 
V sin Is value is clearly incompatible with the spectroscopic de- 
termination of the rotational velocity of the star. 

We also computed the line-bisector behavior following the 
procedure described in Santos et al. (2002) . As shown in Fig. [3] 
a specific signature of the line-bisectors can be found during the 
transit and is anti-correlated with the RVs due to the fact that the 
crossing planet mainly affects the bottom of the spectral lines. 

Assuming the B07 stellar mass and period, of 1.298+;}°^^^0 
and 5.63341 days, and using the fitted ratios Rp/R., and a/Rg, 
we also provide a new determination of the system parameters. 
Those include: orbital semi-major axis, a, stellar radius, Rg, and 
the planetary radius, Rp, mass, Mp, density pp, and surface grav- 
ity gp. Our results are similar to those of B07 with smaller errors. 
The RVs describing the spectroscopic transit are quite noisy and 
some of the measured RVs (around orbital phase 0.885) present 
an unexpected RV shift. This shift does not seem to be due to an 
instrumental deviation and may be explained by a guidance de- 
centring of the telescope. However, the amplitude of the Rossiter 
anomaly is large enough that we can still estimate the parameters 
of this system with quite good uncertainties. 



start Mid End 

0.2 r — I , 1 ^ , 1 — n 1 1 1 • — ; 




Orbital 



Fig. 2. Top: Radial velocities of HAT-P-2 as a function of the 
orbital phase after subtraction of the Keplerian model and super- 
imposed with the best fit of the Rossiter-McLaughlin effect. The 
filled and open circles represent the RV measurements obtained 
during the first and the second sequence of observations of the 
transit with SOPHIE, respectively. The open circle with dotted 
error bars represent one measurement from the Keck set which 
is a few minutes before the ingress phase. Bottom: RV residuals 
after subtracting the orbital solution and the modeled RM effect. 
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Fig. 3. Line-bisectors signature during the spectroscopic transit 
as a function of the orbital phase. As in Fig. |2] the filled and 
open circles represent the RV measurements obtained during the 
first and second sequences of transit observations with SOPHIE, 
respectively. 

5. Discussion and Conclusion 

The value of HAT-P-2b radius (0.95 Rjup) puts this object in 
the mass-radius diagram as an intermediate case between Hot- 
Jupiters and low-mass stars. Its mean density of 12.5 g cm"-' is 
in between the Hot- Jupiter density (0.34 - 1.34 g cm"^) and the 
density of the sr tiallest transiting M dw arfs OGLE-TR-122 and 
OGLE-TR-123 (iPont et al.l l20oi l2()06h. which are 7 5 and 51 g 
cm , respectively. A similar paper (" Winn et al.ll2007h published 
results that are consistent with our conclusions. Super-massive 
planets like HAT-P-2b may constitute a new class of stellar com- 
panions, in between Hot-Jupiters and low-mass stars and near 
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Table 2. Refined parameters of the Keplerian orbital solution 
and the parameters of the HAT-P-2 system, compared to those 
determined in (Bakosetal. 2007, B07). To refers to the peri- 
astron time. The fitted parameters are presented at the top of 
the table. The bottom part of the table lists the system parame- 
ters, derived from the fitted ones. We present in the table both 
the spectroscopic determination and the result from the fit of the 
star's vsinL. 
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Parameter 



Value 



Value from B07 



Period [d] * 
To [HJD] 
e 

w [deg] 
K[ms-'] 
Vo [ms-'] 

a/Rs 
A [deg] 

vsin/j jtM [kms"'] 
vsin/, spictm [kms-' 
i [deg] 
Avkl [ms"'] 
Av5a: [ms"'] 
Av52 [ms"'] 
Stellar jitter [ms"'] 



5.63341 (fixed) 
2454213.4794 +000^3 
5163 +"■'"'2' 
189.92 +[°« 
966.9 + 8.3 
-19855.1 ± 5.8 
0.06891 
10.28 

0.2 !!p 

22.9 tli 
21.3 +1.3 
90 

-328 ± 41 
-19584 ± 17 
-27 ± 12 
17 



[Fe/H] (dex) 

M, [Mo]* 
a [AU] 

Rs [Ro] 
Mp [M,,„p] 

Rp [R/ii;;] 

Pp [g cm"^] 
gp [m s"-] 



Adopted from B07 
' Short term jitter 



0.11 ± 0.10 
-4.75 ± 0.02 

1 9Q8 +0062 

0677 

o.uu; / _0.0017 

1 416 
8.62 J 
0.951 

12 5 +2-^ 

707 +30 

-41 



1039 
1053 



5.63341±0.00013 
2454213.369 ± 0.041 
0.520 + 0.010 
179.3 ± 3.6 
1011 ±38 
N.A. 

0.0684 ± 0.0009 

q 77+1.10 
' ' -0.02 
N.A. 

N.A. 

19.8 ± 1.6 
> 84.6 (20-) 
-380 ± 35 
N.A. 
N.A. 
60 



0.12 ± 
-4.72 : 



0.08 

:0.05 



1 9qo+0.062 

0.0677 ± 0.0014 

1 474+0.042 

9.04 ± 0.50 

0.982!»:?^^ 

11.9!1| 
997+44 
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the Brown dwarf population. HAT-P-2b is the first super-massive 
object around a F8 star for which the exact mass has been deter- 
mined. Such a massive close-in planet is not in agreement with 
the type II migration mechanism, which appears to be more ef- 
ficient for planets around low-mass stars ( Ida & Linll2005l) . This 
could suggest a different formation process for this object, such 
as fragmentation, or interactions between the planet and another 
companion or between the planet and the disk in the evolution 
process. 

A system with an elliptic orbit is expected to move toward 
pseudo-synchronization, with the stellar angular rotation veloc- 
ity tuned to near the an^iilar velocity of the companion at pe- 
riastron passage (IZahnlll977h . Given the Keplerian orbital pa- 
rameters, we computed the angular planetary speed at the peri- 
astron position and found that a stellar radius of 1 .42 Rq implies 
a pseudo-synchronization rotational velocity vsin/, of about 
40 kms"' (Hut 1981, equation 43). Because the observed ro- 
tational velocity of the star is only about 21 kms"', the star is 
definitely not pseudo-synchronized to the planetary orbit. iPeald 
d 19991) indicated that the alignment and the synchronization 
timescales are of the same order of magnitude. Therefore, the 
lack of pseudo-synchronization indicates that the system was 
formed with the stellar spin aligned to the orbital angular mo- 
mentum. 



